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Abstract 

We argue that the tt production asymmetry observed at the tevatron might be sim- 
ply explained if the standard SU{3)c QCD theory is extended to SU{Nc) which is 
spontaneously broken at a scale just above the weak scale. The extended gauge inter- 
actions amplify the radiative QCD contribution to the asymmetry and can potentially 
explain the observations if Nc ~ 5. This explanation suggests a relatively low SU{Nc) 
symmetry breaking scale ~ 0.5 — 1 TeV. We check that such a low SU{Nc) symme- 
try breaking scale is consistent with current collider data. Importantly this scenario 
predicts an abundance of striking phenomena which will be probed at the LHC The 
SU{Nc) model also illustrates the idea that a beyond standard model contribution to 
the ti asymmetry might arise primarily via radiative corrections rather than at tree- 
level. 



1 Introduction 



The CDF collaboration has observed [1] an unexpectedly large forward-backward tt 
production asymmetry at the tevatron. With an integrated luminosity of 5.3 fb~^, the 
asymmetry is found to vary strongly with the invariant mass of the ti pair (M^f). At 
the parton level, the required ti asymmetry in the ti rest frame is observed to be: 

Att = 0.475 ± 0.114 for Ma > 450 GeV 

Att = -0.116 ± 0.153 for Mu < 450 GeV . (1) 

These asymmetries can be compared with the next to leading order QCD prediction 
of[l] [see also [2]]: 

Att = 0.088 ± 0.013 for Mtf > 450 Gey 

Att = 0.040 ± 0.006 for M^- < 450 GeV . (2) 

The statistical significance of the excess is 3.4(7 and represents exciting evidence for new 
physics. Possible explanations for the asymmetry include s-channel interference of an 
exotic color octet gauge boson or t-channel scalar/vector exchange, see e.g. ref.[3, 4] for 
some recent studies. All of these explanations involve an exotic tree-level contribution 
to the ti production amplitude and generally have significant experimental constraints 
coming from resonance production, like-sign top quark production etc. In this paper 
we wish to explore the novel possibility that the asymmetry arises predominately as a 
radiative effect, analogous to the QCD contribution. 

The leading order QCD qq — > ti production asymmetry arises from two differ- 
ent reactions[5]. The first involves interference of initial state with final state gluon 
bremsstrahlung. The second involves interference betweeen the tree-level qq — ?> ti pro- 
cess and its one-loop correction. The Feynman diagrams are shown in Figure 1. 




Figure 1: Interfering qq — > ttj (above) and qq — > tt (below) amplitudes. 
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Considering first the asymmetry due to the interference between the tree-level qq 
tt process and its one- loop correction, it was found in ref.[5] that the QCD asymmetry 
is obtained from the corresponding QED result via the replacement: 

aQEoQ^Qt ^ yJM^""' = 

where Cp = {N^ - l)/2Nc = 4/3 and Tp = 1/2. It is instructive to make the full 
Nc dependence exphcit. Using dlf,^ = {N^ - 4:){N^ - 1)/Nc, we find that the QCD 
asymmetry is obtained from the corresponding QED result via the replacement: 

- 4 

OCQEnQqQt Ols . (4) 

The asymmetry due to the interference of initial state with final state gluon brcmsstrahlung 
is also proportional to d^^^ and has the same Nc dependence. That is, to leading order 
we have that 

- 4 

A,-oc^. (5) 

Thus we arrive at the interesting result that the asymmetry increases rapidly with 
increasing N^.. This observation suggests that gauge models where SU{'i)c is incorpo- 
rated into a larger SU {Nc) gauge symmetry which is broken near the weak scale might 
provide a suitable candidate for the new physics needed to explain the tt asymmetry. 
Furthermore this illustrates the novel idea that a new physics contribution to the tt 
asymmetry could arise primarily via radiative corrections rather than at tree-level. 



2 The SU{b)c model - and generalization to SU{Nc) 



An example of such an extended gauge model for the strong interactions is given by 
the anomaly free SU{5)c ® SU{2)l ® ^/(l)y gauge model proposed some time ago 
by Hernandez and I[6], and further studied in ref.[7, 8, 9, 10, 11]. In the SU{5)c 
model, the quarks and leptons of each generation transform under the gauge symmetry 
SU (5)c <8) SU {2)l ® U{1)y' in the anomaly free representation: 



V 

e 



;i,2,-l), eR~(l,l,-2), 



(6) 



Ql = 



u 

d 



(5, 2, 1/5), UR ~ (5, 1, 6/5), dR - (5, 1, -4/5) . 



(7) 



Gauge symmetry breaking is accomplished via the vacuum expectation values (VEVs) 
of a scalar x ^ (10, 1, 2/5) along with the usual Higgs scalar, 4> ~ (1, 2, -|-1): 
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The SU{5)c ® U{1)y' gauge symmetry is assumed to be spontaneously broken by the 
VEV of X at a scale just above the weak scale, which along with the usual Higgs doublet 
leads to the symmetry breaking pattern: 

SU{5)c^SU{2)l(^U{1)y' 

lix) 

SU{3)c ® SU{2y ® SU{2)l ® U{1)y 

SU{3)c^SU{2y®U{l)Q (9) 

where Y = Y' + is the linear combination of Y' and T5 which annihilates (x) 

(i.e. Y{x) = 0). [T5 is the diagonal generator of SU{5)c orthogonal to those of the 

12", 



SU{3) ® SU{2y subgroup, i.e. T5 = ^fdiag{l/3, 1/3, 1/3, -1/2, -1/2)]. 

The extra colour degree's of freedom have been called 'quirks' and studied in some 
detail by Carlson et a/[ll]. These fermions have electric charge ±1/2 and are confined 
into integer charged bound states by the unbroken and asymptotically free SU (2)' gauge 
interaction [Asu{2y ~ 0.5 GeV]. These exotic fermions are given mass by coupling to 
X via the Yukawa Lagrangian: 

>C = XlQlx{QlT + \RURx{dRY . (10) 

Note that the generation index has been suppressed. Besides these electroweak invari- 
ant masses, there are SU{5)c invariant mass terms coming from the usual coupling 
terms with the standard Higgs doublet, (f). 

In addition to the standard model gauge bosons the model contains a Z'^ and 
charged (with electric charge ±1/6) W'^ gauge bosons (from SU{5) f/(l)y'/S'?7(3) x 
U{1)y coset). The latter transform as a (3, 2, 1/6) and (3, 2, —1/6) under the unbroken 
SU{3)c ® SU{2y ® U{1)q gauge group. These gauge bosons gain masses of[6] 



M^, 5 



5 

(11) 



M|, 12 ' 

where w is the VEV of x- The mass eigenstate neutral 7, Z' gauge bosons couple to 
fermions as follows: 

>C = -eA,J^^ - -^Z,r^ - g,Z'^J%, (12) 
where the currents are given by 

JL = i^Q^i^. Jn = i^N^i;, J^, = i^N'Y^ . (13) 
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The summation of fermion fields is implied and the generators Q, N and N' are given 

by 

Q - is + ^+ 



2 2VT5 ' 
N = (/3-sin2^^g) + A , 

N' = — + A' (14) 
where A and A' are small corrections given by 

48^>2cos2^^' 2VT5glcos^ey, ' ^ ' 

Note that we have defined 6w by cos^ 6^, = m^/m^. 

Importantly the scale of SU (5)c ® U (l)y' symmetry breaking can be very low[6, 7, 
8, 9, 10, 11]. The Z' is very weakly constrained because it couples mainly to hadrons 
and its mixing with the Z is naturally very small. For example, in ref.[7] only the 
very modest limit of mz' ~ 100 GeV was obtained from neutral current phenomenol- 
ogy. A more stringent limit of Mz' > 280 GeV was obtained from UA2 dijet data 
in ref.[9]. In section 4 we will confront this model with more recent dijet data from 
the Tevatron where we show that a low SU{5)c symmetry breaking scale is still ex- 
perimentally allowed, although not without some tension. The quirk mass is also very 
weakly constrained. A mass limit on the quirks can be obtained from measurements of 
the Z boson width, which limits the mass of the quirks to be greater than 43 GeV. In 
addition, flavour changing neutral current processes, such as the radiative contribution 
to K'^ — K'^ mass mixing constrain the quirk masses to be roughly degenerate or the 
corresponding CKM-type matrix to be almost diagonal[ll]. 

The case of larger (odd) Nc is a straightforward generalization to the SU{5)c case, 
but requires a more complex Higgs sector to achieve the required symmetry breaking^. 
To break the symmetry in a phenomenologically consistent way in SU (iVg) with > 5 
will require several scalar multiplets. For example, an SU{l)c®SU{2)L®U{l)yf model 
can be constructed by placing the fermions into the anomaly-free representation: 



= ] ~ (1,2,-1), eK~(l,l,-2), (16) 




Ql ^ \ ^] -(7,2,1/7), ii«~(7,l,8/7), ci«^(5,l,-6/7) . (17) 
The SU{7)c®U{1)yi can be spontaneously broken by the VEVs of two scalars, X11X2 ~ 



^Thc case of even Nc is possible, although the simplest such theories run into phenomenological 
problems [12]. 
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(21, 1, 2/7) at a scale just above the weak scale. These VEVs are given by: 
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(18) 



These VEVs will give the exotic quirk degree's of freedom electroweak invariant masses, 
analogously to the SU{5)c case [Eq.(lO)]. These VEVs along with (0) lead to the 
symmetry breaking pattern: 



SU{7),0SUi2)L®U{l)Y' 
i (Xi), (X2) 

SU{3)c ® SU{2y ® SU{2y ® SU{2)l ® U{1)y 

;(</>) 

SU{3)c (8) SU{2y (8) SU{2y (8) U{1)q 
where Y = Y' + -f^T-j. Here Ty denotes the SU{7) diagonal generator: 



(19) 



24 

= J-diag{l/3, 1/3, 1/3, -1/4, -1/4, -1/4, -1/4) . 



(20) 



The SU (7)c model will contain two Z' gauge bosons, Z[, Z'^. These gauge bosons couple 
to the generators 



z' 



nl , . /I/ Tx 

cosO — + sme* — 
2 2 

T T 
— sm fc* — + cos — 



(21) 



where T7 is defined above and = \/2(iia5f(0, 0, 0, — 1/2, — 1/2, 1/2, 1/2), and 9' is 
obtained by diagonalizing the Z' gauge boson mass matrix: 



12^2 



f^sK + ^i) \/y92[wi - ^2 J 



That is, tan 2^' = 4^ 



. For SU{7)c and the case of larger Nc, other symmetry 



breaking patterns are possible depending on the scalar content of the theory. 
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3 Top quark forward-backward asymmetry in SU (Nc) 
color theories 



The QCD contributions to the tt forward-backward asymmetry feature infrared (IR) 
divergences. The contribution due to gluon bremsstrahlung can be separated into 

hard and soft gluon emission, and the IR divergence from the latter cancelling the 
IR divergence from the box diagam contribution to the asymmetry [2, 5]. Calculations 
show that the hard gluon bremsstrahlung contribution to Af^ is negative, while the soft 
gluon plus virtual part is positive [2, 5]. It turns out that the soft gluon bremsstrahlung 
+ virtual box diagram part is about twice the magnitude of the hard bremstrahlung 
part, so that the predicted QCD asymmetry is positive. 

In the extended SU {Nc) color gauge models, the tt asymmetry should be approx- 
imately the same as that predicted in the standard model at energies below the scale 
of SU {Nc) symmetry breaking, which we take here for illustration as a single scale, m. 
The radiative processes involving exotic W, Z' gauge bosons and quirks are suppressed 
by powers of VS/m. Thus, things are clear in the m ^ limit - the standard model 
calculation of the tt asymmetry results. Let us now examine the m <^ limit. In this 
limit we expect the SU{Nc) symmetry to be a good approximation and provided that 
T quirks (we denote the quirk partner of the corresponding quark with an uppercase 
letter) decay into t-quarks ^, then the tt asymmetery should follow from Eq.(5). If the 
region M^f > 450 GeV is sufficiently high that the SU (Nc) limit becomes useful, then 
we can predict from Eq.(5) a tt forward-backward asymmetry of: 

Att ^ 0.22 for Nc^5 

Att ^ 0.34forA^e = 7 

Att ~ 0.45 for A^c = 9 

Att ~ 0.56 for Nc = 11 etc. (22) 

In the realistic case we might have \/^ ~ m and that T quirks might not decay pre- 
dominately into t-quarks, in any case the bremsstrahlung contribution can potentially 
be kinematically supressed. In this case, the exotic contribution to the tt asymmetry 
should arise predominately from the interference of the tree-level qq tt amplitude 
with the virtual box diagram amplitude. It seems that the size of this contribution 
can be even larger than that suggested by the SU{Nc) symmetry estimate, Eq.(22), 
since that estimate includes the hard bremstrahlung contribution, which is negative in 
sign. It is also possible that the exotic box diagram contributions are smaller than the 
SU{Nc) estimate, Eq.(22) or even that it has the wrong sign (although it is unlikely 
that it has the wrong sign for all ranges of parameters). For completeness one also 
needs to calculate the box diagram contribution with gauge bosons replaced by the x 
scalar /Goldstone bosons (in 't Hooft-Feynman gauge). The latter could dominate if the 

^The part of the asymmetry due to the interference of initial and final state bremsstrahlung will 
contain qq tW'T final states (as well as qq — )• tig final states). The T quirks can decay into top 
quarks: T tUu via a virtual W , assuming that the U quirk is the lightest quirk. 
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quirks are heavier than the gauge bosons. In this paper we will be content with point- 
ing out the novel possibility that the top quark forward-backward asymmetry might be 
radiatively induced with SU (Nc) models as a promising example. Obviously if the CDF 
anomaly is confirmed by future data then more detailed calculations will be necessary 
to pin down the mquirk/f^z',w' l^c parameter space allowed by the experiments. 

Observe that the total top quark production cross section in SU {N^) color models is 
not expected to be significantly 15%) modified from the standard model prediction. 
This is because the main effect of the new physics is to increase the radiative contribu- 
tion which generates the asymmetry via interference with the QCD contribution. The 
proportion of T quirks produced is relatively small. 

If SU (Nc) color gauge theory is the origin of the tt asymmetry then the scale of new 
physics is likely to be low (unless N^. is very big), and we thus expect this explanation 
to be tested in the future at the LHC experiment. However in the meantime we need to 
check that the low SU (Nc) symmetry breaking scale is compatible with existing data. 



4 Tevatron limits on SU (5)c [SU {Nc)] symmetry break- 
ing scale 



Data from the Tevatron might constrain the SU{5)c symmetry breaking scale. The 
most promising signature is expected to be the effects of the Z' since this particle 
can be produced on resonance if its mass is low[8, 9]. The Z' will manifest itself 
primarily as a bump in the dijet invariant mass distribution. No such bump is seen, 
leading to strong constraints on possible Z' couphng to quarks. For example, the CDF 
collaboration [13] are able to exclude an axigluon with mass 260 < m < 1250 GeV. 
However, it turns out that the SU{5)c Z' couples to each quark color with coupling 
gs/\/T5- The l/\/T5 factor significantly suppresses the cross section c.f. axigluons. We 
have computed the cross section for pp — > Z' — > 2jets. This is obtained by calculating 
the cross section for the parton subprocess qq ^ Z' ^ q'q' and then folding this in 
with the parton distribution functions in the usual way. We use the 2008MSTW NLO 
parton distribution functions in our numerical work[14]. The cross section is multiplied 
by the standard K factors 



to incorporate the QCD corrections in the initial and final states. We have also applied 
a rapidity cut of \y\ < 1 for each of the two jets. Our results are shown in figure 2. 
The calculation assumes that the quirks are light, so that they contribute to the total 
width of the Z' but do not contribute to the dijet signal'^. Also shown in the figure is 
the corresponding CDF 95% C.L. upper hmit[13]. 

^Quirk pairs can annihilate into jets, but at a lower invariant mass, where the dijet constraints are 
much weaker. 



Snasiq'^) 
9 




TT 



(23) 



7 




400 600 800 1000 1200 1400 

Z' Mass [GeV] 

Figure 2: The Z' particle production cross section times branching fraction into dijets 
times the acceptence for both jets to have \y\ < 1. Also shown, is the corresponding 
CDF 95% C.L. upper limit. 



As the figure shows, there is some tension between a low SU{5)c breaking scale, 
which disfavours the model but is not quite strong enough to conclusively exclude it. 
In the case of SU{7)c and larger color groups, having multiple Z's might help evade 
the CDF dijet constraint by effectively smoothing out the dijet bump. The Z' can 
also decay leptonically via a small Y' coupling [from the A' part in Eqs.(14, 15)]. 
However the branching fraction to leptons is very small, for example we find that 
Br{Z' — )■ jijl) ~ 1.4 X 10~^. In fact this branching fraction is sufficiently small to evade 
tevatron limits[15] on resonances decaying into muon pairs. 

The tevatron might potentially constrain the quirk masses. Quirk pairs can be 
produced via 7, Z, Z' exchange. While heavy quirks can decay into lighter ones, the 
lightest pair of quirks arc expected to hadronize (assuming the lightest quirk is lighter 
than the W or SU{2y colored scalar components). As discussed in some detail in 
ref . [1 1] , this quirk pair can be viewed as a highly excited quirkonium state. Deexcitation 
occurs via numerous emission of SU{2)' glueballs, dubed 'hueballs' in ref. [11]. Note 
that the hueball emission will lead to substantial missing energy but only small missing 
transverse momentum. The quirkonium will end up in the ground state where it will 
decay via virtual 7, Z into standard model particles. 

The quirk production cross section times branching ratio into muons at the tevatron 
has been calculated in ref. [11]. They obtain a cross section times branching ratio into 
muon pairs of about 10/6 for a quirk mass of 100 GeV for left-handed quirks, and much 
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weaker limits for right-handed quirks. This is well below current limits obtained by 
the CDF collaboration [15]. Thus, we conclude that a low SU{Nc) symmetry breaking 
scale of a few hundred GeV is phenomenologically viable, and consequently the SU{Nc) 
explanation of the top quark forward-backward asymmetry is possible. 

Observe that the SU{Nc) explanation of the tt forward-backward asymmetry would 
imply a similarly sized asymmetry for bb quark production. It has been argued in 
ref.[4, 16] that such an asymmetry, although challenging, could hkely be measured at 
the tevatron. If this is the case then this will further test the SU (Nc) model. 

The CDF collaboration has also recently found an anomalously high number of 
boosted jets (pt > 400 GeV for the leading jet) with invariant mass near the top quark 
mass[17]. The anomaly could be due to an underestimation of background but might 
also be a signal for new physics. Such a signal can be explicable in SU (N^) models. 
For instance, if mz' ~ 0.8 — 1 TeV, then Z' — > tt decays will generally lead to a boosted 
top. Furthermore, from figure 2 we see that the cross section for hadronically decaying 
Z' is of order 0.1 pb, which implies a cross section to hadronically decaying top pairs 
of about 10 fb. This value is roughly consistent with the value necessary to explain the 
CDF data if it is due to new physics[17, 18]. The CDF collaboration did not observe 
any anomaly in the boosted jet plus missing energy channel (corresponding to one of 
the top quarks decaying semileptonically) . However it has been argued in ref.[19] that 
the semileptonic channel is less sensitive than the fully hadronic decay channel. 

The LHC should be able to probe the SU{Nc) explanation for the tt asymmetry. 
The quirks, W, Z' gauge bosons and also the exotic colored x scalar can be produced 
at the LHC, leading to dijet signals. We expect, however, that at least several jb~^ of 
data might need to be accumulated before this explanation of the tt asymmetry can 
be put to the test. We leave detailed predictions for the LHC for future work. 

5 Conclusion 

In conclusion, we have pointed out that the tt forward-backward asymmetry observed at 
the tevatron might be simply explained if the standard SU (3)c QCD theory is extended 
to SU (Nf.) which is spontaneously broken at a scale just above the weak scale. The 
extended gauge interactions amplify the QCD radiative contribution to the asymmetry 
and can potentially explain the observations if Nc 5. This explanation suggests a 
relatively low SU{Nc) symmetry breaking scale ^ 0.5 — 1 TeV, and we have checked 
that this is compatible with recent collider data. Importantly this scenario predicts 
an abundance of new physics just above the weak scale, which will be probed in the 
future at the LHC. 

Acknowledgements 

The author thanks Archil Kobakhidze for useful conversations and to G. Rodrigo for 
some helpful correspondence. This work was supported in part by the Australian 
Research Council. 



9 



References 

[1] T. Aaltonen et al, (CDF collaboration), arXiv: 1101.0034 (2011). 

[2] L. G. Almeida, G. F. Sterman and W. Vogelsang, Phys. Rev. D78: 014008 (2008); 
O. Antuano, J. H. Kuhn and G. V. Rodrigo, Phys. Rev. D77: 014003 (2008); 
M. T. Bowen, S. D. Ellis and D. Rainwater, Phys. Rev. D73: 014008 (2006); For 
earlier work see also R. W. Brown, D. Sadhev and K. O. Mikaelian, Phys. Rev. 
Lett. 43: 1069 (1979); F. Halzen, P. Hoyer and C. S. Kim, Phys. Lett. B195: 74 
(1987). 

[3] P. Ferrario and G. Rodrigo, Phys. Rev. D78: 094018 (2008) [arXiv: 0809.3354]. S. 
Jung et al, Phys. Rev. D81: 015004 (2010) [arXiv: 0907.4112]; K. Cheung, W-Y. 
Keung and T-C. Yuan, Phys. Lett. B682, 287 (2009) [arXiv: 0908.2589]; A. Arhrib, 
R. Benbrik, C-H. Chen, Phys. Rev. D82: 030434 (2010) [arXiv: 0911.4875]; P. H. 
Frampton, J. Shu and K. Wang, Phys. Lett. B683, 294 (2010) [arXiv: 0911.2955]; 
D-W. Jung et al, Phys. Lett. B691, 238 (2010) [arXiv: 0912.1105]; J. Cao et 
al, Phys. Rev. D81: 014016 (2010) [arXiv: 0912.1447]; I. Dorsnet et al, Phys. 
Rev. D81: 05509 (2010) [arXiv: 0912.097] J. Shu, T. M. P. Tait and K. Wang, 
Phys. Rev. D81: 034012 [arXiv: 0911.3237]; V. Barger, W-Y. Kueng, C-T. Yu, 
Phys. Rev. D81: 113009 (2010) [arXiv: 1002.1048]; arXiv: 1102.0279; Q-H. Cao 
et al, Phys. Rev. D81: 114004 (2010) [arXiv: 1003.3461]; B. Xiao, Y. Wang and 
S. Zhu, Phys. Rcv.D82: 034026 (2010) [arXiv: 1006.2510]; C-H. Chen, G. Cvetic 
and C. S. Kim, Phys. Lett. B694, 393 (2011) [arXiv: 1009.4165]; K. Cheung, T- 
C. Yuan, arXiv: 1101.1445; E. L. Berger et al, arXiv: 1101.5625; J. Cao et 
al, arXiv :1101.4456; C. Delaunay et al, arXiv: 1101.2902; B. Grinstein, A. L. 
Kagan, M. Trott and J. Zupan, arXiv: 1102.3374; K. M. Patel and P. Sharma, 
arXiv: 1102.4736; B. Bhattacherjee, S. S. Biswal and D. Ghosh, arXiv: 1102.0545; 
E. Ramirez Barreto, Y. A. Coutinho and J. Sa Borges, arXiv: 1103.1266; A. R. 
Zerwekh, arXiv: 1103.0956. 

[4] Y. Bai, J. L. Hewett, J. Kaplan and T. G. Rizzo, arXiv: 1101.5203. 

[5] J. H. Kuhn and G. Rodrigo, Phys. Rev. D59: 054017 (1999) [arXiv: hep- 
ph/9807420]. 

[6] R. Foot and O. F. Hernandez, Phys.Rev.D41:2283 (1990), Erratum-ibid.D42:948 
(1990). 

[7] S. L. Glashow and U. Sarid, Phys.Lett.B246:188 (1990). 
[8] R. Foot, O. F. Hernandez and T. G. Rizzo, Phys.Lett.B246:183 (1990). 
[9] R. Foot, O. F. Hernandez and T. G. Rizzo, Phys.Lett.B261:153 (1991). 
[10] O. Hernandez, Phys. Rev. D44: 1997 (1991). 



10 



[11] E. D. Carlson, L. J. Hall, U. Sarid, J. W. Burton, Pliys.Rev.D44: 1555 (1991). 

[12] R. Foot and O. F. Hernandez, Mod. Phys. Lett. A7: 1315 (1992). 

[13] T. Aaltonen et al (CDF Collaboration), Phys. Rev. D79: 112002 (2009) [arXiv: 
0812.4036]. 

[14] A. D. Martin, W. J. Stirling, R. S. Thorne, G. Watt, Eur. Phys. J. C63: 189 
(2009) [arXiv: 0901.0002]. 

[15] T. Aaltonen et al (CDF Collaboration), arXiv: 1101. 4578. 

[16] M. J. Strassler, arXiv: 1102.0736. 

[17] CDF Collaboration, CDF pubhc note 10234. 

[18] K. Blum et al, arXiv: 1102.3133. 

[19] Y. Eshel, O. Gedaha, G. Perez and Y. Soreq, arXiv: 1101.2898. 



11 



